Dust Properties in the FUV in Ophiuchus 
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ABSTRACT 

We have derived the albedo (a) and phase function asymmetry factor (g) 
of interstellar dust grains at 1100 A using archival Voyager observations of 
diffuse radiation in Ophiuchus. We have found that the grains are highly forward 
scattering with g = 0.55 ± 0.25 and a = 0.40 ± 0.10. Even though most of the gas 
in this direction is in the Ophiuchus molecular cloud, the diffuse FUV radiation 
is almost entirely due to scattering in a relatively thin foreground cloud. This 
suggests that one cannot assume that the UV background is directly correlated 
with the total amount of gas in any direction. 

Subject headings: dust — ultraviolet: ISM — infrared: ISM 



1. INTRODUCTION 

The Ophiuchus molecular cloud is one 
of a large complex of clouds at a distance 
of about 160 pc from the Sun (Chen et al. 
1997). The region has been extensively 
studied in CO (de Geus et al. 1990; Dame 
et al. 2001) as well as with four-colour pho- 
tometry (Corradi et al. 2004) allowing us 
to determine the three-dimensional distri- 
bution of the matter in this direction. Our 



interest in this region is that one of the 
first observations of diffuse emission in the 
far ultraviolet (FUV - 912 - 1216 A) was 
made here by Holberg (1990) who iden- 
tified the emission as starlight from the 
nearby Scorpius-Centaurus OB association 
scattered by dust in the Ophiuchus molec- 
ular cloud. 

Because of the technical difficulties in- 
herent in diffuse observations in the FUV, 
largely due to scattering from the intense 
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geocoronal Lyct line, there have been very 
few observations of the background radi- 
ation in this spectral region (see Bowyer 
1991; Henry 1991, for reviews of the dif- 
fuse radiation in both the near and far ul- 
traviolet). By far the largest and most re- 
liable data set has come from the ultravi- 
olet spectrographs (UVS) aboard the two 
Voyager spacecraft which made observa- 
tions of various astrophysical targets dur- 
ing the interplanetary phase of their mis- 
sion, between their hugely successful plan- 
etary encounters. Many of these targets 
were of objects with no intrinsic FUV flux 
and Murthy et al. (1999) used them for a 
comprehensive study of the diffuse FUV 
radiation field. Because of the sensitiv- 
ity of the Voyager UVS to diffuse radiation 
and the distance of the spacecraft far from 
the Earth where emission from interplan- 
etary H I is minimized, these remain the 
definitive observations of the diffuse radia- 
tion field in the FUV. 

Except for a small extragalactic com- 
ponent at high latitudes (Henry 1991), 
the diffuse UV radiation is largely due to 
starlight scattered by interstellar dust and 
as such can be used to derive the scattering 
properties of the interstellar dust grains. 
Unfortunately there has been considerable 
controversy about both the level of the dif- 
fuse radiation and the modeling used to 
extract the optical constants (see Draine 
2003a; Mathis et al. 2002, for a discus- 
sion of the difficulties) and there have been 
only loose observational constraints on the 
albedo (a) and phase function asymmetry 
factor (g) of the dust grains. 

As mentioned above, Holberg (1990) 
discovered intense diffuse emission from 
Ophiuchus using the Voyager UVS and in- 



terpreted this emission as starlight back- 
scattered by the Ophiuchus molecular 
cloud. We have reexamined these obser- 
vations along with others from the Voy- 
ager archives (Murthy et al. 1999) in the 
light of a more sophisticated model and an 
improved understanding of the dust distri- 
bution in the direction of Ophiuchus (Cor- 
radi et al. 2004) and have found that the 
emission is actually due to the forward 
scattering of the light from a much lower 
density sheet of material in front of the 
Ophiuchus molecular cloud. There were 
many observations spread throughout the 
entire region and so we have been able, 
for the first time, to remove the degener- 
acy between the albedo (a) and the phase 
function asymmetry parameter (g) which 
have plagued studies of the diffuse radia- 
tion (Draine 2003a). We have constrained 
o to 0.40 ± 0.10 and g to 0.55 ± 0.25, in 
good agreement with the theoretical pre- 
diction of Weingartner & Draine (2001) for 
a mixture of graphitic and silicate grains. 

2. OBSERVATIONS 

The two Voyager spacecraft include 
identical Wadsworth-mounted objective 
grating ultraviolet spectrographs (UVS) 
covering the wavelength region between 
500 and 1700 A, with their greatest sen- 
sitivity at wavelengths below 1200 A. 
The field of view of the spectrographs is 
0.1° x 0.87° with a spectral resolution of 
38 A for diffuse sources. The spacecraft 
were launched within two weeks of each 
other - - Voyager 2 in 1977 August and 
Voyager 1 in 1977 September — and ob- 
tained a wealth of information on all four 
of the giant planets. The two spacecraft 
are still continuing operation at the edge 
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of the solar system with more than 10,000 
days of operation each. A full description 
of the UVS spectrographs and further in- 
formation about the interstellar mission of 
the Voyager spacecraft is given by Holberg 
& Watkins (1992). 

While the spacecraft were between plan- 
etary encounters, they observed many as- 
tronomical targets. Amongst these tar- 
gets were a series of scans in the vicinity 
of Ophiuchus with the Voyager 2 space- 
craft in 1982 (Holberg 1990). We have sup- 
plemented these with further observations 
from the Voyager archives taken at vari- 
ous times between 1982 and 1994 (Murthy 
et al. 1999). There were a total of 31 such 
locations and these are plotted on an IRAS 
100 fim map in Fig. 1 and tabulated in Ta- 
ble 1 with the total hydrogen column den- 
sity (Schlegel et al. 1998) and the IRAS 
100 jum flux (Wheelock et al. 1994). There 
was a limit cycle motion of a few tenths of 
a degree in the Voyager pointing and the 
position reported in the table is an average 
of the actual pointing of the instrument. 

The data used here are from Murthy 
et al. (1999) and are available from the 
authors. The data reduction is fully de- 
scribed in that paper and in Holberg (1986) 
and consists of fitting three components 
to the observed signal: dark noise from 
the spacecraft's radioisotope thermoelec- 
tric generator (RTG); emission lines from 
interplanetary H I; and emission from cos- 
mic sources. The dark noise was sub- 
tracted using the continuum below the Ly- 
man limit 912 A and the astronomical and 
heliospheric emission were then simultane- 
ously fit using templates for the emission. 
This reduction procedure was shown to be 
consistent and reproducible over observa- 



tions separated widely in time and between 
the two spacecraft and yielded la limits 
as low as 30 photons cm~ 2 s _1 sr _1 A -1 
over a large part of the sky. For the bright 
sources reported on here, the diffuse emis- 
sion dominated the raw data and there was 
little uncertainty in the derived levels. 

3. MODEL 

The diffuse emission in Ophiuchus is the 
result of light from the nearby stars scat- 
tering from interstellar dust in the line of 
sight and were modeled in a similar manner 
to Shalima & Murthy (2004) in the Coal- 
sack. The stellar radiation field was calcu- 
lated at the location of the scattering dust 
using the model of Sujatha et al. (2004), 
in which the distance and spectral type of 
each star was taken from Hipparcos data 
(Perryman et al. 1997) and the flux calcu- 
lated using Kurucz models (Kurucz 1992) 
with the latest modifications from his web 
page (http://kurucz.harvard.edu). Eight 
stars (see Fig. 1) contribute ~90% of the 
total ISRF in this region. The properties 
of the stars are given in Table 2. 

The amount of radiation scattered to 
the observer is dependent on the scattering 
function of the grains and we have used the 
Henyey-Greenstein scattering phase func- 
tion, cj)(9) (Henyey & Greenstein 1941). 



YK 1 Att (1 + g 2 - 2gcos(6)) 1 - 5 v 1 

Here a is the albedo — which can range 
from for dark grains to 1 for fully reflect- 
ing grains — and g is the phase function 
asymmetry factor with g = indicating 
isotropic scattering and higher values indi- 
cating forward scattering grains. 
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Draine (2003b) has suggested that the 
Henyey-Greenstein function underestimates 
the scattered radiation for highly forward 
scattering grains (g > 0.7) in the FUV. We 
have found, empirically that using the the- 
oretical scattering function of Weingartner 
& Draine (2001) for a mixture of graphite 
and silicate grains makes no more than a 
10% difference in the derivation of the opti- 
cal constants for g = 0.65 and so, for con- 
sistency with the literature, we only cite 
the results using the Henyey-Greenstein 
function. The primary uncertainty in our 
procedure is in the location of the scatter- 
ing dust and we discuss this below. 

The interstellar medium in the direction 
of Ophiuchus is dominated by the huge 
molecular complexes in Ophiuchus (Fig. 1) 
at a distance of ~ 160 pc (Chen et al. 
1997). Because of the thickness of the 
cloud, only foreground stars will contribute 
to the diffuse radiation observed at the 
Earth. Back-scattering from the molecular 
cloud is an order of magnitude too small to 
account for the observed UV intensity and 
hence we conclude that the scattering must 
be from two extended sheet-like structures 
which Corradi et al. (2004) have shown to 
cover the entire region between the galactic 
longitudes of 290° - 10° and latitudes of 
-25° — 1-25°. The nearer of the two sheets 
is at a distance of d < 60 pc from the Sun 
and the other, from which most of our ob- 
served scattering comes, lies between 100 
and 150 pc from the Sun, depending on the 
direction. The latter sheet is likely part 
of the neutral ring surrounding the com- 
plex of molecular clouds in Ophiuchus dis- 
covered by Egger & Aschenbach (1995) in 
the ROSAT all-sky survey. Corradi et al. 
(2004) has found a column density of 3.2 



x 10 19 cm -2 in the nearer cloud and a 
much larger column density of 3.7 - 27 
x 10 20 cm -2 in the further cloud. Other 
than these two clouds and the Ophiuchus 
molecular cloud, there is very little mate- 
rial out to a distance of at least 200 pc from 
the Sun (Frisch et al. 2005). 

In our model, we have divided the to- 
tal H I from Dickey & Lockman (1990) in 
any line of sight into the two foreground 
clouds, with a constant value of 3.2 x 10 19 
cm~ 2 in the nearer cloud and the rest in the 
further cloud. For the Ophiuchus molecu- 
lar clouds, we have subtracted the H I col- 
umn density from the total hydrogen col- 
umn density of Schlegel et al. (1998) and 
distributed this excess material at the loca- 
tion of the molecular cloud. We converted 
the H I column density to a dust scatter- 
ing cross-section using the theoretical val- 
ues of Weingartner & Draine (2001), im- 
plicit in which is the dust to gas ratio of 
Bohlin et al. (1978). In practice, the ob- 
served UV emission is almost entirely from 
the more distant of the two sheets and so is 
most sensitive to the exact distance of that 
sheet, or rather to the distance between the 
sheet and the stars dominating the ISRF in 
this region, and the amount of dust in that 
cloud. Because this distance is uncertain, 
we have used a 3 parameter model in which 
we allow the distance to the dust to vary 
but fix a and g to a common value over all 
31 positions. We then use a single scatter- 
ing model to calculate the scattered flux 
in the UV and independently calculate the 
thermal emission at 100 fim for the best 
fit parameters. The best fit distances are 
plotted in Fig. 2 and show a variation of 
100 - 125 pc in this region, consistent with 
the values found by Corradi et al. (2004). 
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4. RESULTS AND DISCUSSION 

It is interesting to note that the level of 
the diffuse emission is not at all correlated 
with the amount of material in the line of 
sight (Fig. 3). Instead there is a tight cor- 
relation between the level of the ISRF and 
the scattered radiation (Fig. 4). This has 
important implications for the study of the 
diffuse radiation field. It is often claimed 
that the diffuse radiation is correlated with 
the H I column density (e.g. Bowyer 1991; 
Schiminovich et al. 2001); however, it is 
becoming increasingly evident that the dif- 
fuse radiation is dominated by local effects, 
particularly near bright stars (c.f. Murthy 
& Sahnow 2004). 

Our model predictions match the obser- 
vations extremely well both in UV scatter- 
ing (Fig. 5) and in 100 /xm infrared emis- 
sion (Fig. 6). We have derived a 90% con- 
fidence contour (as per Lampton, Margon, 
& Bowyer 1976) for a and g and this is 
shown in Fig. 7. Our 90% confidence lim- 
its on a and g are 0.40 ± 0.10 and 0.55 ± 
0.25 respectively and are consistent with 
the theoretical predictions of Weingartner 
& Draine (2001) for average Milky Way 
dust with Ry=3.1. 

There are very few determinations of 
the optical parameters of interstellar grains 
at wavelengths shorter than 1200 A (Draine 
2003a; Gordon 2004) and most of these 
have come from observations of reflec- 
tion nebulae. Witt et al. (1993) found 
an albedo of 0.42 ± 0.04 from Voyager 2 
observations of NGC 7023 and Burgh et al. 
(2002) found an albedo of 0.30 ± 0.10 from 
rocket observations of NGC 2023. Both 
groups claimed that the grains were highly 
forward scattering with g ~ 0.8. Shalima 



& Murthy (2004) found a similar value of 
0.40 ± 0.20 through Voyager observations 
of the Coalsack Nebula but were not able 
to constrain g. 

Although our derived albedo (0.40 ± 
0.10) is consistent with the earlier deter- 
minations, we find a slightly lower value 
for g (0.55 ± 0.25). It is possible that con- 
ditions are different in reflection nebulae 
as opposed to the diffuse ISM we observe 
or it may be that, as Draine (2003b) sug- 
gests, the scattering function is poorly rep- 
resented by the Henyey-Greenstein func- 
tion in the FUV leading to differences be- 
tween determinations in different geome- 
tries, particularly for highly forward scat- 
tering grains. 

5. CONCLUSIONS 

We have used Voyager observations of 
the diffuse FUV radiation in the region of 
Ophiuchus to investigate the optical con- 
stants of the interstellar dust grains. We 
have found that the intense emission in 
this region arises not in the dense molecu- 
lar cloud which contains most of the mat- 
ter in this direction but rather in a much 
thinner neutral sheet in front of the cloud. 
In fact, there is no correlation between ei- 
ther the 100 /im IRAS emission or N(H) 
and the scattered FUV light. Instead the 
FUV emission is tightly correlated with 
the strength of the local ISRF. Thus un- 
like thermal dust emission in the IR, where 
the dust is optically thin and the material 
along the entire line of sight contributes to 
the total, scattering in the FUV requires 
much thinner clouds and nearby bright 
stars. 

In general, one should not expect cor- 
relations on global scales between the UV 
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and IR or UV and N(H). This is borne out 
by Murthy & Sahnow (2004) who find only 
a mild dependence on the 100 /im for the 
FUV flux but not by Schiminovich et al. 
(2001) who claim a correlation between the 
UV (at 1500 A) and N(H) using data from 
the NU VIEWS rocket flight. We are pur- 
suing further investigations with data from 
the Galaxy Evolution Explorer (GALEX) 
to test these correlations. 

We find that the interstellar dust grains 
are highly forward scattering with a g (— 
{cos9)) of 0.55 ± 0.25 and an albedo of 0.40 
± 0.10. These results are in general agree- 
ment both with the theoretical predictions 
of Weingartner & Draine (2001) and previ- 
ous observations of reflection nebulae. As 
a test of our model, we have independently 
calculated the 100 ^m intensities for each 
of the locations for our best fit a and g 
values and found them to match the IRAS 
(100 fim) observations. 
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Fig. 1. — IRAS 100 jum map of the region 
is shown with contours labelled in units of 
MJy sr _1 . The filled squares show the loca- 
tions of the Voyager observations and the 
asterisks show the positions of the bright- 
est UV stars in the region. 
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Fig. 3. — Observed UV intensity at each 
location is plotted against the correspond- 
ing values of the total, N(H) (dotted line) 
and neutral, N(H I) (solid line) hydrogen 
column densities. There is clearly no cor- 
relation between the FUV intensity and ei- 
ther N(H) or N(H I). 
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Fig. 2. — Best fit distances to the scat- 
tering layer for each of the locations as a 
function of galactic longitude in degrees. 
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Fig. 4. — Observed UV intensity at each 
location is plotted against the ISRF. The 
non-zero intercept is due to the absorption 
of the diffuse radiation in the intervening 
ISM. 
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Fig. 6. — Modeled IR intensities corre- 
sponding to an albedo of 0.40 and g = 
0.6 have been plotted against the observed 
IRAS (100 jum) values for each location. 
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Fig. 5.— Modeled FUV (1100 A) inten- 
sities corresponding to a = 0.40 and g = 
0.6 have been plotted against the observed 
values for each location. 
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Fig. 7. — 90% confidence contour for all 
31 locations (solid contour) is shown. The 
contour corresponds to a limit of 0.40 ± 
0.10 and 0.55 ± 0.25 on the albedo and 
g respectively. Also plotted is the inter- 
section of the individual 90% confidence 
contours for each of the locations (dashed 
line) . 
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Table 1 

Details of observed locations in Ophiuchus 



Location 


1 


b 


Observed FluxiError 


IRAS(100 fim) a 


N(H) 6 




(deg) 


(deg) 


(ph cm~ 2 s _1 sr _1 A -1 ) 


(MJy sr" 1 ) 


(x 10 21 cm- 2 ) 


1 


359.8 


17.8 


1900±170 


43 


3.3 


2 


0.2 


22.1 


2200±140 


38 


3.8 


3 


7.9 


7.8 


2710±440 


21 


2.5 


4 


7.8 


8.4 


3090±260 


17 


2.4 


5 


0.6 


21.4 


3200±140 


36 


2.4 


6 


1.4 


19.7 


3270±360 


30 


2.8 


7 


0.4 


16.8 


3320±140 


37 


3.8 


8 


1.8 


16.2 


3320±180 


32 


3.1 


9 


355.5 


31 


3510±50 


19 


1.1 


10 


5.8 


21.6 


3510±75 


33 


2.5 


11 


2.4 


19.2 


3740±310 


26 


2.1 


12 


348.1 


12.1 


4080±360 


40 


2.1 


13 


2.4 


18.7 


4310±65 


32 


2.8 


14 


4.1 


16.8 


4430±295 


30 


4.7 


15 


355.5 


30.3 


4620±90 


16 


1.1 


16 


2.4 


26.7 


4710±150 


27 


1.8 


17 


3.6 


17 


4990±50 


28 


4.5 


18 


3.2 


18.2 


5100±90 


30 


3.1 


19 


5.9 


14 


5120±360 


28 


2.6 


20 


350.4 


9.8 


5580±325 


42 


2.4 


21 


348.5 


12.2 


5760±445 


40 


2.0 


22 


359.2 


23.5 


5940±95 


33 


2.4 


23 


357.5 


27.5 


6520±125 


18 


1.2 


24 


356.4 


29.2 


7020±260 


15 


1.2 


25 


358.1 


25.2 


7340±110 


26 


1.4 


26 


357 


28.3 


8840±380 


18 


1.1 


27 


10.2 


22.4 


8920±610 


25 


3.3 


28 


9.1 


23.1 


9700±680 


27 


2.7 


29 


9.7 


22.5 


10250±540 


25 


2.6 


30 


354.3 


13 


10900±570 


49 


2.6 


31 


350.5 


24.9 


33000±955 


25 


1.4 



a Wheclock et al. (1994). 
b Column densities from Schlegel et al. (1998). 



10 



Table 2 

Properties of contributing stars 



Star 


1 


b 


Spectral Type 


Temperature 


log (g) 


log (z) 




(deg) 


(deg) 




(K) 






COph 


6.28 


23.59 


09V 


35000 


3.94 





X Oph 


357.93 


20.68 


B2Vne 


22000 


3.94 





v Sco A 


354.61 


22.70 


B2IV 


22000 


3.94 





u> Sco 


352.75 


22.77 


B1V 


25600 


3.94 





8 Sco 


350.10 


22.49 


B0.2IV 


30000 


3.94 





IT SCO 


347.22 


20.23 


B1V 


25600 


3.94 





r Sco 


351.54 


12.81 


B0V 


30000 


3.94 





A Lib 


350.72 


25.38 


B3V 


19000 


3.94 
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